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Abstract In this proceedings we present a state-of-the-art method of calculating 
thermodynamic potential at finite temperature and finite chemical potential, using 
Hard Thermal Loop perturbation theory (HTLpt) up to next-to-next-leading-order 
(NNLO). The resulting thermodynamic potential enables us to evaluate different 
thermodynamic quantities including pressure and various quark number suscepti¬ 
bilities (QNS). Comparison between our analytic results for those thermodynamic 
quantities with the available lattice data shows a good agreement. 


1 Introduction 

With the increasing advancements of relativistic heavy ion colliders, the study of 
the color deconfined quark gluon plasma (QGP) phase has also become more ex¬ 
tensive. From the last couple of decades, paired with the recurring successes of 
non-perturbative methods like lattice QCD (1QCD) simulations, one unavoidable 
question also haunted the theoreticians, now and again. How far can we trust those 
ideas which describes the QGP, based on perturbation theory? To answer this ques¬ 
tion, along with the high loop order calculations now one is also interested in the 
comparison of those analytic results with the relevant lattice QCD data. 
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QCD pressure in perturbation theory for both zero and finite chemical potentials 
has already been calculated and results are known up to [1-3]. The re¬ 

sults show a poor convergence in the phenomenologically important region near the 
transition temperature. Search for the root of this problem eventually suggests some 
systematic reorganization of bare perturbation theory. One of them, in which we 
will focus on, is the Hard Thermal Loop perturbation theory (HTLpt), developed by 
Andersen, Braaten and Strickland [4,5] using the concept of Hard Thermal Loops, 
proposed by Braaten and Pisarski [6], almost two decades ago. Since then, HTLpt 
has become one of the primary tools for evaluating thermodynamic quantities and 
it has already been used in the evaluation of QCD thermodynamic potentials up to 
one-loop (LO) [7-12], two-loop (NLO) [13-15] and recently even up to three-loop 
(NNLO) [16-19] for both zero and finite chemical potentials [20,21]. 

In this proceedings, we present the recent evaluation of the NNLO QCD ther¬ 
modynamic quantities [19,20] in HTLpt at finite temperature and finite chemical 
potentials both for flavor-independent and flavor-dependent cases. We emphasize 
that this NNLO calculation of thermodynamic potential is completely analytic and 
gauge independent. 


2 Results 


In the present section we show some plots showing the behavior of the thermody¬ 
namic quantities based on the final results from [19,20]. Though we used both the 
one-loop and three-loop running coupling, but in this proceedings we show results 
only with the former one. By taking the value of a s (1.5GeV) = 0.326 from lattice 
measurements [22], the fixed value of the scale comes out to be 176 MeV. 
For purely gluonic and fermionic Feynman diagrams, two separate renormalization 
scales, A g and A, are used respectively. We use the debye mass prescription intro¬ 
duced by Braaten and Nieto in [23] which is generalized to finite chemical potential 
by Vuorinen in [2,3]. 

In Fig. 1 we show the scaled NNLO HTLpt pressure for both Hb = 0 and iLn = 400 
MeV, which can be obtained directly from the NNLO thermodynamic potential. 
In these plots along with the others to follow, the thick black line represents the 
result obtained using the central values of renormalization scales, e.g. A g = 2nT 
and A = 2 K \/T 2 + fx 2 / K 2 . The surrounding light-blue band specifies the variation 
of the result by varying these scales with a factor of two, e.g. nT <A g < 4nT. The 
NNLO pressure is compared with available 1QCD data which shows a very good 
agreement [24,26]. 

Usually the baryon number susceptibility (BNS) is related to the quark number 
susceptibilities [19,20]. If one assumes flavor independent unanimous quark chem¬ 
ical potential (ji u = j. 1 ,/ = )X S = [i ), then BNS can be calculated in a straightforward 
way by taking the derivatives of the pressure with respect to /i. In Fig.2 we com¬ 
pare the scaled second-order (left) and fourth-order (right) BNS with various lattice 
data [25,27,28,30,31]. 
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Fig. 1 Scaled NNLO HTLpt pressure for Nf = 2 + 1 at )l B = 0 (left) and fig = 400 MeV (right), 
compared with lattice data from Borsnyi et al. [24,26] 



T [MeV] T [MeV| 


Fig. 2 For Nf — 2 + 1, scaled NNLO HTLpt second order (left) and fourth order (right) baryon 
number susceptibility compared with various lattice data, labeled as WB, BNL-BI(B). BNL- 
BI(u,s), MILC, and TIFR come from Refs. [25], [27], [28], [30], and [31], respectively. 


To compute single quark number susceptibilities (sQNS), flavor dependent quark 
chemical potentials are taken into account in the general expression of NNLO ther¬ 
modynamic potential [20]. Derivatives with respect to /!,- yield either diagonal (same 
flavor on all indices) or off-diagonal (different flavor on some or all indices) suscep¬ 
tibilities. In Fig.3 we compare NNLO HTLpt fourth-order diagonal sQNS (left) and 
the only non-vanishing fourth-order off-diagonal sQNS (right) with available 1QCD 
data from various groups [27-29,31]. 


3 Conclusions 

In this proceedings, the NNLO QCD thermodynamic potential and various thermo¬ 
dynamic quantities within 3-loop HTLpt at finite temperature and chemical poten- 
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Fig. 3 For Nf = 2 + 1, NNLO HTLpt fourth order diagonal single quark number susceptibility 
(left) and the only non-vanishing fourth order off-diagonal quark number susceptibility (right), 
compared with various lattice data, labeled as BNL-BI(uudd), BNL-BI(u, s), BNL-BI(uuss), and 
TIFR taken from Refs. [27], [28], [29], and [31], respectively. In the left figure the dashed blue line 
indicates the Stefan-Boltzmann limit for this quantity. 


tial(s) are presented. Results obtained using the central values of the renormalization 
scales are well in accord with the accessible lattice data, for some of the quantities 
even down to phenomenologically more important temperatures relevant to LHC. 
On the ground of these perturbative results, one may further study other important 
facets of Quark Gluon Plasma, e.g. transport properties. 
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